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(S| ■ ABSTRACT 

It is shown how deviations of neutrino flavor mix from the canonical 1:1:1 
«— h ' in high energy astrophysical neutrinos can provide evidence for neutrino decays 

^5 , with lifetimes as long as 10 4 s/eV. 

cr\ : 

1. Introduction 

In my talk at the last Venice meeting, I discussed neutrino fluxes, and especially 
flavor composition of the neutrino beams, from astrophysical sources of high energy 
cn ' neutrinos. It is generally agreed that pion and muon decays produced in beam dump 

type situations are the most likely sources of neutrinos in these environments^. The 
cn ■ fraction of heavy flavors, which can produce u' T s, even at these high energies, is very 

smalf^. In this case the initial flavor mix is given to a good approximation by 
&, e :0„„:^«l:2:O. 

With the current knowledge of 5m?-, and for the large L/E that we anticipate, the 
oscillations average out so that on arrival at earth, the flux of each mass eigenstate 
are given by <j> Vi = Y, a 4>t°J lTce I U a i | 2 , where U a i are elements of the neutrino mixing 
.^J | matrix. For three active neutrino species (as we assume throughout) there is now 

strong evidence to suggest that u u and u T are maximally mixed and ~ 0. The 
consequent — v T symmetry means that in the mass eigenstate basis the neutrinos 
are produced in the ratios <f> Vl : 4> U2 : <p U3 = 1:1:1, independent of the solar mixing 
angle. Oscillations do not change these proportions, but only the relative phases 
between mass eigenstates, which will be lost. An incoherent mixture in the ratios 1 
: 1 : 1 in the mass basis implies an equal mixtu re in any basis, and in particular the 
flavor basis in which the neutrinos are detected I^M). 

There is also the possibility that the initial flavor mix is not the canonical 1:2:0. 
For example, the muons can lose energy (as at high energies in the atmosphere) e.g. 
in a magnetic field; in that case no high energy v' e s are available and the flavor mix 
becomes ^ (j) Ve : (p Ufi : <p UT = 0:1:0. The effect of mixing then is to produce a 
final mixture which is given by (fi Ue : (p u ^ : 4> Ut = a : 1 : 1 with a between 1/2 to 2/3 
depending on the precise value of the solar mixing angle. 

This shows the value of measuring the flavor mix; it can yield information on 
the production mechanism and the environment; even in the presence of mixing and 



oscillations. But it can do even more; such as yield information on neutrino properties. 
There are other ways in which the flavor mix can change from the canonical 1:1:1 
figure. One in particular, which gives rise to striking signatures, is the decay of 
neutrinos^. Before discussing these, let me motivate the possibility of neutrino decay. 

2. Neutrino Decaj^ 

We now know that neutrinos have non-zero masses and non-trivial mixings. This 
is based primarily on the evidence for neutrino mixings and oscillations from the data 
on atmospheric neutrinos and on solar neutrinos. 

If this is true then in general, the heavier neutrinos are expected to decay into 
the lighter ones via flavor changing processes. The only questions are (a) whether 
the lifetimes are short enough to be phenomenologically interesting (or are they too 
long?) and (b) what are the dominant decay modes. 

Throughout the following discussion, to be specific, I will assume that the neutrino 
masses are at most of order of a few eV. Since we are interested in decay modes which 
are likely to have rates (or lead to lifetimes) which are phenomenologically interesting, 
we can rule out several classes of decay nodes. 

First, consider radiative decays, such as v e — > Vj + 7. Since the experimental 
bounds on fi Ui , the magnetic moments of neutrinos, come from reactions such as 
v e e — > e'V" which are not sensitive to the final state neutrinos; the bounds apply to 
both diagonal as well as transition magnetic moments and so can be used to limit 
the corresponding lifetimes. The bounds should really be on mass eigenstateJ^, but 
since the mixing angles are large, it does not matter much. The current bounds ar 

r Ue > 5.10 18 sec 
r Vfl > 5.10 16 sec 

t Vt > 2.10 11 sec. (1) 

There is one caveat in deducing these bounds. Namely, the form factors are 
evaluated at q 2 ~ 0(eV 2 ) in the decay matrix elements whereas in the scattering 
from which the bounds are derived, they are evaluated at q 2 ~ 0(MeV 2 ). Thus, some 
extrapolation is necessary. It can be argued that, barring some bizzare behaviour, 
this is justifiecPl. 

An invisible decay mode with no new particles is the three body decay vi — > UjUji>j. 
Even at the full strength of Z coupling, this yields a lifetime of 2.10 34 s, far too long 
to be of interest. There is an indirect bound from Z decays which is weaker but still 
yields 2.10 30 sElJ. 

Thus, the only decay modes which can have interestingly fast decays rates are two 
body modes such as z/j — > Vj + x and z/j — ► Vj + x where x is a very light or massless 
particle, e.g. a Majoron. 



The only possibility for fast invisible decays of neutrinos seems to lie with Majoron 
models'' 7 '. There are two classes of models; the 1=1 Gelmini-RoncadeHf^ majoron 
and the 1=0 Chikasige-Mohapatra-Peccei*'*"' majoron. In general, one can choose the 
majoron to be a mixture of the two; furthermore the coupling can be to flavor as well 
as sterile neutrinos. The effective interaction is of the form: 



g a V0V a J (2) 

giving rise to decay: 

u a -> vp (or vp) + J (3) 

where J is a massless J = L = 2 particle; v a and v$ are mass eigenstates which 
may be mixtures of flavor and sterile neutrinos. Models of this kind which can giv e 
rise to fast neutrino decays and satisfy the bounds below have been discussed^. 
These models are unconstrained by fi and r decays which do not arise due to the 
AL = 2 nature of the coupling. The 1=1 coupling is constrained by the bound on the 
invisible Z width; and requires that the Majoron be a mixture of 1=1 and 1= (P. 
The couplings of and v e (g^ and g e ) are constrained by the limits on multi-body 
7T, K decays tc — > {ivvv and K — > {ivvv and on /i — e university violation in 7r and 
K decays*" 2 ', but not sufficiently strongly to rule out fast decays. Direct limits on 
such decay modes are also very weak. The current strongest bound come from solar 
neutrino results and are of order oir/m > lCr 4 s/ This limit is based primarily 

on the non-distortion of the Super-Kamiokande spectrurrP-^ , and takes into account 
the potentially competing distortion caused by oscillations (see also Ref!*^) as well 
as the appearance of active daughter neutrinos. The SN 1987A data place no limit 
at all on these neutrino decay modes, since decay of the lightest mass eigenstate is 
kinematically forbidden, and with the inherent uncertanties in the neutrino fluxes, 
even a reasonable V\ flux alone can account for the datcP^. 

The strongest lifetime limit is thus too weak to eliminate the possibility of astro- 
physical neutrino decay by a factor about 10 7 x (L/100 Mpc) x(10 TeV/E). Some 
aspects of the decay of high-energy astrophysical neutrinos have been considered in 
the past. It has been noted that the disappearance of all states except v\ would 
prepare a beam that could in principle be used to measure elements of the neutrino 
mixing matrix, namely the ratios : : U^^ ] . The possibility of measuring 

neutrino lifetimes over long baselines was mentioned in Ref.E2J, and some predictions 
for decay in four-neutrino models were given in Ref!^. We have shown that the 
particular values and small uncertainties on the neutrino mixing parameters allow 
for the first time very distinctive signatures of the effects of neutrino decay on the 
detected flavor ratios. The expected increase in neutrino lifetime sensitivity (and cor- 
responding anomalous neutrino couplings) by several orders of magnitude makes for a 
very interesting test of physics beyond the Standard Model; a discovery would mean 
physics much more exotic than neutrino mass and mixing alone. As shown below, 



neutrino decay because of its unique signature cannot be mimicked by either different 
neutrino flavor ratios at the source or other non-standard neutrino interactions. 

A characteristic feature of decay is its strong energy dependence: exp^—Lm/Er), 
where r is the rest-frame lifetime. For simplicity, we will assume that decays are 
always complete, i.e., that these exponential factors vanish. The assumption of com- 
plete decay means we do not have to consider the distance and intensity distributions 
of sources. We assume an isotropic diffuse flux of high-energy astrophysical neutrinos, 
and can thus neglect the angular deflection of daughter neutrinos from the trajectories 
of their parents^' . 

Disappearance only. Consider the case of no detectable decay products, that 
is, the neutrinos simply disappear. This limit is interesting for decay to 'invisible' 
daughters, such as a sterile neutrino, and also for decay to active daughters if the 
source spectrum falls sufficiently steeply with energy. In the latter case, the flux of 
daughters of degraded energy will make a negligible contribution to the total flux at 
a given energy. Since coherence will be lost we have. 

<f>» a = E €°; rce (E) | Up | 2 | U ai | 2 e- L '^ (4) 

if} 

l>5 £ ^^{E)\U pi \ 2 \U ai \\ (5) 

instable), f3 

where the <p Ua are the fluxes of u a , U a i are elements of the neutrino mixing matrix 
and t are the neutrino lifetimes in the laboratory frame. Eq. (5) corresponds to the 
case where decay is complete by the time the neutrinos reach us, so only the stable 
states are included in the sum. 

The simplest case (and the most generic expectation) is a normal hierarchy in 
which both z/3 and z/2 decay, lea ving only the lightest stable eigenstate v\. In this case 
the flavor ratio is : : U*P l K Thus if U e3 = 

<Pue ■ <t> v » ■ 4>u T = cos 2 6 Q : - sin 2 6 & : - sin 2 6 & ~ 6 : 1 : 1, (6) 

where 8 & is the solar neutrino mixing angle, which we have set to 30°. Note that this 
is an extreme deviation of the flavor ratio from that in the absence of decays. It is 
difficult to imagine other mechanisms that would lead to such a high ratio of v e to 
v^. In the case of inverted hierarchy, z/3 is the lightest and hence stable state, and so 

K ■ K ■ K = U 2 e3 ■■ U% :( % 0:1:1. (7) 

If U e3 = and 6 atm = 45°, each mass eigenstate has equal and v T components. 
Therefore, decay cannot break the equality between the 0j, M and 4> Ut fluxes and thus 
the <fi Ue : <f> v ratio contains all the useful information. The effect of a non-zero U e3 
on the no-decay case of 1 : 1 : 1 is negligible. For the normal hierarchy, the decay 



flavor mix can change z/ e /V M from 6/1 to 3/1 and u^/ur from 1/1 to 3/1. In the case 
of inverted hierarchy, there is no such dependence. 

Appearance of daughter neutrinos. - If neutrino masses are quasi-degenerate, 
the daughter neutrino carries nearly the full energy of the parent. An interesting and 
convenient feature of this case is that we can treat the effects of the daughters without 
making any assumptions about the source spectra. Including daughters of full energy, 
we have 

ME) ^ E <f>f urce ( E ) I Uf* l 2 l u <* I 2 ( 8 ) 

i(3 

+E0r rce (£)i f/ «i 2 i f/ » i 2 

where B is a branching fraction and stable and unstable states are denoted henceforth 
by % and j respectively. 

If instead the neutrino mass spectrum is hierarchical, the daughter neutrinos will 
be degraded in energy with respect to the parent, so that 

K(E) ^? E €T rce (E) I Ufa | 2 | U m | 2 (9) 

i/3 

dE> W E , E J2€7 rce (E') I U M | 2 | U ai | 2 B^, 

HP 

where E is the daughter and E' is the parent energy. The normalized energy spectrum 
of the daughter is given by 

1 dT(E'.E) . . 

W f >f = (10) 

EE F(E') dE { } 

If the neutrinos are Majorana particles, daughters of both helicities will be detectable 
(as neutrinos or antineutrinos), whereas if they are Dirac particles, daughters of one 
helicity will be sterile and hence undetectable. In the hierarchical limit, the energy 
distributions in the lab frame are E/E" 2 and (£" — E) / E" 2 respectively for the helicity 
conserving and helicity flipping modes. 

In the case of Majorana neutrinos, we may drop the distinction between neutrino 
and antineutrino daughters and sum over helicities. Assuming the source spectrum 
to be a simple power law, E~ a , we find 

4> Va {E) L ^ E €°; rce (E) I Upi | 2 | U m | 2 (11) 

i(3 

+ W€7 Ce (E)\U^\ 2 \U m | 2 E,_n 



This is identical to the expression in Eq. (8) except for the overall factor of 1/a in front 
of the second term. For Dirac neutrinos we detect only the daughters that conserve 



helicity, the effect of which is only to change the numerical coefficient of the second 
sum in Eq. (11). Thus, although the flavor ratio will differ from the cases above, it 
is still independent of energy-i.e., decay does not introduce a spectral distortion of 
the power law. We stress that we have assumed a simple but reasonable power law 
spectrum E~ a ; a broken power law spectrum, e.g. would lead to a more complicated 
energy dependence. 

Uniqueness of decay signatures.-Depending on which of the mass eigenstates 
are unstable, the decay branching ratios, and the hierarchy of the neutrino mass 
eigenstates, quite different ratios result. For the normal hierarchy, some possibilities 
are shown in Tabled 



Table 1: Flavor ratios for various decay scenarios. 
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The most natural possibility with unstable neutrinos is that the heaviest two mass 
eigenstates both completely decay. The resulting flavor ratio is just that of the lightest 
mass eigenstate, independent of energy and whether daughters are detected or not. 
For normal and inverted hierarchies we obtained 6:1:1 and 0:1:1 respectively. 
Interestingly, both cases have extremes Ve : (fi u ^ ratios, which provides a very useful 
diagnostic. Assuming no new physics besides decay, a ratio greater that I suggests 
the normal hierarchy, while a ratio smaller than 1 suggests an inverted hierarchy. In 
the case that decays are not complete these trends still hold. 

An important issue is how unique decay signatures would be. Are there other 
scenarios (either non-standard astrophysics or neutrino properties) that would give 
similar ratios? There exist astrophysical neutrino production models with different 
initial flavor ratios, such as : 1 : ^ , for which the detected flavor ratios (in the 
absence of decay) would be about 0.5 : 1 : 1. However, since the mixing angles 9 Q 
and d a tra are both large, and since the neutrinos are produced and detected in flavor 
states, no initial flavor ratio can result in a measured <p Ve : <fi v ratio anything like 
that of the two main cases, 6:1:1 and : 1 : 1. 

In terms of non-standard particle physics, decay is unique in the sense that it 
is "one way", unlike, say, oscillations or magnetic moment transitions. Since the 
initial flux ratio in the mass basis is 1 : 1 : 1, magnetic moment transitions between 



(Majorna) mass eigenstates cannot alter this ratio, due to the symmetry between 
i — > j and j — > i transitions. On the other hand, if neutrinos have Dirac masses, 
magnetic moment transitions (both diagonal and off-diagonal) turn active neutrinos 
into sterile states, so the same symmetry is not present. However, the process will 
not be complete in the same way as decay-it will average out at 1/2, so there is no 
way would be left with a only single mass eigenstate. 

Experimental Detectability. Deviations of the flavor ratios from 1:1:1 due 
to possible decays are so extreme that they should be readily identifiable^. Up- 
coming high energy neutrino experiments, such as Ice cubd^, will not have perfect 
abilities to separarely measure the neutrino flux in each flavor. However, the quan- 
tities we need are closely related to observables, in particular in the limit of — v T 
symmetry (8 atm = 45° and U e3 = 0), in which all mass eigenstates contain equal 
fractions of and v T . In that limit, the fluxes for and v T are always in the ratio 1 
: 1, with or without decay. This is useful since the v T flux is the hardest to measure. 

Detectors such as IceCube will be able to directly measure the flux by long- 
ranging muons which leave tracks through the detector. The charged-current interac- 
tions of v e produce electromagnetic showers. However, these are be hard to distinguish 
from hadronic showers caused by all flowers through their neutral-current interactions, 
or from the charged-current interactions of v T (an initial hadronic shower followed by 
either an electromagnetic or hadronic shower from the tau lepton decay We thus 
consider our only experimental information to be the number of muon tracks and the 
number of showers. 

The relative number of shower events to track events can be related to the most 
interesting quantity for testing decay scenarios, i.e., the v e to ratio. The precision 
of the upcoming experiments should be good enough to test such extreme flavor ratios 
produced by decays. If electromagnetic and hadronic showers can be separated, then 
the precision will be even better. 

Comparing, for example, the standard flavor ratios of 1 : 1 : 1 to the possible 
6:1:1 generated by decay, the more numerous electron neutrino flux will result 
in a substantial increase in the relative number of shower events. The details of 
this observation depends on the range of muons generated in or around the detector 
and the ratio of charged to neutral current cross sections. This measurement will 
be limited by the energy resolution of the detector and the ability to reduce the 
atmospheric neutrino background. The atmospheric background drops rapidly with 
energy and should be negligibly small at and above the PeV scale. 

Discussion and Conclusions.- The flux ratios we discuss are energy-independent 
because we have assumed that the ratios at production are energy- independent, that 
all oscillations are averaged out, and that all possible decays are complete. The 
first two assumptions are rather generic, and the third is a reasonable simplifying 
assumption. In the standard scenario with only oscillations, the final flux ratios are 
(p Ue : (j) v : 4>v T = 1:1:1- In the cases with decay, we have found rather different 



possible flux ratios, for example 6 : 1 : 1 in the normal hierarchy and : 1 : 1 in the 
inverted hierarchy. These deviations from 1:1:1 are so extreme that they should 
be readily measurable. 

If we are very fortunatd^, we may be able to observe a reasonable number of 
events from several sources (of known distance) and/or over a sufficient range in 
energy. Then the resulting dependence of the flux ratio iy e /vp) on L/E as it evolves 
from say 6 (or 0) to 1, can be clear evidence of decay and further can pin down the 
actual lifetime instead of just placing a bound. 

To summarize, we suggest that if future measurements of the flavor mix at earth 
of high energy astrophysical neutrinos find it to be 

K/Kl^r = a/1/1; (12) 

then 

(i) a = 1 (the most boring case) confirms our knowledge of the MN!5p^ matrix and 

our prejudice about production mechanism; 

(ii) a ~ 1/2 indicates that the source emits pure v'^s and the mixing is conventional; 

(iii) a > 1 indicates that neutrinos are decaying with normal hierarchy; and 

(iv) o; <C 1 would mean that neutrino decays are occuring with inverted hierarchy. 
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